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ABSTRACT 

We use absolutely calibrated data from the ARCADE 2 flight in July 2006 to model Galactic 
emission at frequencies 3, 8, and 10 GHz. The spatial structure in the data is consistent with a 
superposition of free-free and synchrotron emission. Emission with spatial morphology traced by the 
Haslam 408 MHz survey has spectral index /3synch = — 2.5 ± 0.1, with free- free emission contributing 
0.10 ± 0.01 of the total Galactic plane emission in the lowest ARCADE 2 band at 3.15 GHz. We 
estimate the total Galactic emission toward the polar caps using either a simple plane-parallel model 
with csc(|6|) dependence or a model of high-latitude radio emission traced by the COBE/FIRAS map 
of Cli emission. Both methods are consistent with a single power-law over the frequency range 22 
MHz to 10 GHz, with total Galactic emission towards the north polar cap Tcai — 0.498 ± 0.028 K and 
spectral index f} = —2.55 ± 0.03 at reference frequency 1 GHz. The well calibrated ARCADE 2 maps 
provide a new test for spinning dust emission, based on the integrated intensity of emission from the 
Galactic plane instead of cross-correlations with the thermal dust spatial morphology. The Galactic 
plane intensity measured by ARCADE 2 is fainter than predicted by models without spinning dust, 
and is consistent with spinning dust contributing 0.4 ± 0.1 of the Galactic plane emission at 22 GHz. 
Subject headings: radio continuum: ISM, radiation mechanisms: non-thermal, cosmic microwave 
background 



1. INTRODUCTION 

The cosmic microwave background (CMB) is a valu- 
able probe of physical conditions in the early universe. 
Its frequency spectrum records the history of energy 
transfer between the evolving matter and radiation fields 
to constrain the energetics of the early universe. We view 
the CMB through diffuse emission from the interstel- 
lar medium. At centimeter wavelengths, the dominant 
contributions are from synchrotron emission originat- 
ing from electrons accelerated in the Galactic magnetic 
field, and free-free emission (thermal bremsstrahlung) 
from electron-ion collisions. Thermal emission from in- 
terstellar dust is negligible, but electric dipole emis- 
sion from a population of small, rapidly rotating dust 
grains could contribute a substantial fraction of the 
total Galactic emission at wavelengths near 1 cm 
■ jDraine & Lazarian' '1998'; 'Dob ler fc E^inkbeiner I l2007l : 
' iMiv illc-Deschenes et_al. 2008). Absolute measurements 
of the sky temperature at centimeter wavelengths mea- 
sure diffuse emission to separate Galactic from primor- 
dial emission and provide information on physical pro- 
cesses in the interstellar medium. 

The Absolute Radiometer for Cosmology, Astro- 
physics, and Diffuse Emission (ARCADE) is an instru- 
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ment to measure the absolute temperature of the sky in 
search of distortions from a blackbody spectrum. AR- 
CADE operates at centimeter wavelengths between full- 
sky surveys at radio frequencies below 3 GHz and the Far 
Infrared Absolute Spectrophotometer (FIRAS) survey at 
frequencies above 60 GHz. It consists of a set of cryogenic 
radiometers observing at 37 km altitude from a balloon 
pay load. Each radiometer uses a double- nulled design, 
measuring the temperature difference between a corru- 
gated horn antenna and an internal reference as the an- 
tenna alternately views the sky or a full-aperture black- 
body calibrator. The internal reference can be adjusted 
to null the antenna-reference signal difference, while the 
calibrator temperature can be independently adjusted 
to null the sky-calibrator signal difference. ARCADE 
thus measures small shifts about a precise blackbody, 
greatly reducing dependence on instrument calibration 
and stability. The calibrator, antennas, internal refer- 
ence, and radiometer front-end amplifiers are mounted 
within a large liquid helium Dewar and are maintained 
near thermal equilibrium with the CMB. Boiloff helium 
vapor, vented through the aperture, forms a barrier be- 
tween the instrument and the atmosphere to allow op- 
eration in full cryogenic mode with no windows between 
the optics and the sky. 

A 2-channel prototype (ARCADE 1) fiew in 2001 
and again in 2003, observing the sky at 10 and 30 
GHz to demon strate the feasibility of cryogenic open - 
aperture optics (jKogut et al. ll2004l : [Fixsen et al. Il2004( l. 
A second-generation 6-channel instrument (ARCADE 2) 
flew in 2005 and again in 2006. A motor fail ure in 2005 al- 
lowed only a single channel to view the sky (jSingal et al. I 
l2006f ). The 2006 flight successfully obtained observa- 
tions at 3.3, 8.3, 10.2, 30, and 90 GHz. Data from the 
2006 flight are consistent with a blackbody CMB spec- 
trum, but show a clear detection of an extragalactic ra- 
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Fig. 1. — Sky maps from the ARCADE 2 2006 flight (thermodynamic temperature in Galactic coordinates), 
masked in grey. The Galactic plane is clearly visible. 
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dio background. This paper describes the analysis of 
Galactic emission from the 2006 flight. C ompanion pa- 
pers describe the ARCADE 2 instrument (jSingal et al. I 
[2008i) ■ the calibrat ion and sky temperature analysis 
(jFixsen et al. II2008I ) and the impl ications of the detect ed 
extra-Galactic radio background (jSeiffert et al. |[200^ . 

2. OBSERVATIONS AND SKY MAPS 

The ARCADE 2 instrument launched from the 
Columbia Scientific Balloon Facility in Palestine, TX car- 
rying a complement of 7 Dicke-switched radiometers and 
1800 liters of liquid helium. During ascent, the Dicke 
switch failed in one radiometer. The remaining radiome- 
ters began sky observations on 2006 July 22 at 05:08 UT, 
continuing through 08:11 UT just prior to flight termi- 
nation. 

Figure [1] shows the absolute sky temperatures binned 
by Galactic coordinates. Each radiometer views the sky 
through a corrugated conical feed horn, scaled with wave- 
length to produce a Gaussian beam w ith 11? 6 full width 
at half maximum (jSingal et al. l[2005l ). The beams point 
30° from the zenith and the entire payload spins at 0.6 
RPM so that the beams scan a circle of 60° diameter cen- 
tered on the zenith. The spin frequency is comparable 
to the 1// knee of the instrument noise, so that strip- 
ing at the mK level is visible in the binned data. Each 
radiometer observed roughly 7% of the sky. 



The measured sky temperatures are dominated by the 
CMB monopole. Galactic emission is clearly visible in 
the data at 3, 8, and 10 GHz. At higher frequencies, 
the smaller Galactic signal and higher instrument noise 
combine to prevent a clear detection of Galactic emission. 
Throughout this paper, we consider only data from the 3, 
8, and 10 GHz radiometers. Table [T] summarizes the in- 
flight performance of each radiometer. We compute the 
instrument noise by comparing the variance of the data 
within each sky bin to the number of observations in that 
bin. This value for the noise thus includes calibration as 
well as the effect of low-frequency striping. 

The observed emission is a combination of CMB, syn- 
chrotron, and free-free emission. Synchrotron emission 
results from the acceleration of cosmic ray electrons in 



TABLE 1 

Frequency Bands For Foreground Analysis 



Frequency 
(GHz) 


Bandwidth 
(MHz) 


White Noise 
mK sl/2 


3.15 


210 


11.8 


3.41 


220 


10.1 


7.98 


350 


5.5 


8.33 


350 


5.2 


9.72 


860 


3.0 


10.49 


860 


3.0 
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the Galactic magnetic field. For a power-law distribu- 
tion of electron energies N{E) oc E^p propagating in a 
uniform magnetic field, the synchrotron emission is also 
a power law, 

TAil^) (X (1) 

where Ta is antenna temperatur e, v is the radiation 
frequ ency, and /3s — ~{p + 3)/2 (jRvbicki fc LightmaiTI 
Il979| ). Free-free emission is also a power law in antenna 
temperature, with spectral index /3s = —2.15. 

We evaluate the observed spectral index in the AR- 
CADE 2 sky maps by plotting the temperature observed 
at one frequency against the temperature of the same sky 
pixel observed at a different frequency. Since the AR- 
CADE 2 absolute calibration is derived from comparison 
of sky data to a blackbody calibrator of known phys- 
ical temperature, the sky s olutions have units of ther- 
modynamic temperature (iFixsen et a l. 2008). In order 
to evaluate the spectral index, we convert the sky maps 
(Fig. [T]) from thermodynamic to antenna temperature 

Ta = T (2) 

exp(a;) — 1 

where T is thermodynamic temperature, x = hv/kT, 
h is Planck's constant, and k is Boltzmann's constant. 
Figure [2] shows the resulting TT plots. 

The same Galactic features are observed at each fre- 
quency. A linear fit to each plot yields spectral index 
P = -2.43 ± 0.03 from 3.3 to 8.3 GHz, ^ = -2.47 ± 0.02 
from 3.3 to 10.3 GHz, and f3 = -3.12 ± 0.17 from 8.3 
to 10.3 GHz, where the quoted errors include the mea- 
surement uncertainties in both maps for each plot. The 
measured index is consistent with a superposition of syn- 
chrotron and free-free emission. 

3. SPATIAL STRUCTURE 

A search for distortions in the CMB spectrum requires 
correction of Galactic foreground emission. The AR- 
CADE 2 data has only six frequency channels contain- 
ing significant Galactic emission (Table [1]), and can not 
uniquely separate the observed emission into individual 
components (synchrotron, free- free, or other sources). 
We combine the ARCADE 2 data with other sky sur- 
veys in order to model foreground emission. 

A widely-used technique describes the observed emis- 
sion at each frequency i/ as a linear combination of fixed 
"template" maps, 

r^(^.,p)=^a,(^.)G,(p) (3) 

i 

where p is a pixel index and Gi (p) are a set of sky maps 
tracing different components of the interstellar medium. 
The spectral dependence of the fitted coefficients a(i^) 
can then be ascribed to the component traced by each 
template. 

Several choices for temp late maps t r acing free-free 
emission are available. iGold et al. I (|2008 ^ use a 
maximum-entropy algorithm to derive a map of free-free 
emission based on the Wilkinson Microwave Anisotropy 
Probe (WM AP) 5-yea r data and a compilation of Ha 
emission (Fi nkbeinerl l2003f l. Since the maximum- 
entropy method (MEM) template is heavily weighted by 
the microwave data in the Galactic plane, it minimizes 
extinction effects in the Ha map. As alternatives to the 
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Fig. 2.— TT plots from the ARCADE 2 sky maps in units of 
antenna temperature. The best-fit spectral index (solid line) is con- 
sistent with a superposition of synchrotron and free-free emission. 

maximum-entropy templ ate, we also use either a map 
of thermal dust emission ("Finkbe iner. Davis, fc Schlegel I 
([l999 ) model 8 ) or the FIRAS m ap of C ii emission at 158 
fxm (jFixsen. Bennett, fc Mather1ll999f ). neither of which 
is affected by extinction. 

The choice of synchrotron template depends on the 
treatment of spatial variation in the synchrotron spectral 
index. To the extent th at this variation can b e neglected, 
the 408 MHz survey (jHaslam et al. I Il981f ) provides a 
high signal-to-noise ratio tracer for synchrotron emission. 
Since the 408 MHz survey contains both free-free and 
synchrotron emission, we correct it for free-free emission 
using the MEM free-free template scaled to 408 MHz us- 
ing a spectral index f3s = —2.15 so that to first order 
the template represents only synchrotron emission. We 
convolve both the synchrotron and free-free templates to 
match the ARCADE 11?6 beam width. 

The ARCADE 2 data are dominated by the CMB 
monopole, even towards the Galactic plane. Both the 
Galactic and extragalactic signals also contain monopole 
terms, as do the synchrotron and free-free templates. We 
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TABLE 2 

Coefficients for Spatial Templates 



tVequency 


Off 


Ols 


(GHz) 


(mK nW"-"^ sr) 


(mK/K) 


3.15 


3.22 ± 0.11 


2.02 ± 0.05 


3.41 


3.06 ± 0.09 


1.70 ± 0.04 


7.98 


0.26 ± 0.05 


0.24 dz 0.02 


8.33 


0.28 ± 0.05 


0.24 dz 0.02 


9.72 


0.39 ± 0.03 


0.04 ± 0.01 


10.49 


0.37 ± 0.03 


0.05 ± 0.01 



thus include a monopole template as a nuisance parame- 
ter and compute the template coefficients aiiy) by min- 
imizing 

2 _ \p [Ta{p, f) - aoji^) - as(iy)Gs{p) - as{i')Gs{p)f 

(4) 

where ao, ag, and as are the coefficients for the 
monopole, free-free, and synchrotron templates, respec- 
tively, and cr(p, v) is the instrument noise in each pixel 
and frequency channel. 

The spatial structure in the ARCADE 2 data may be 
described as a superposition of emission traced by the 
Haslam 408 MHz survey and the FIRAS map of Cii emis- 
sion. Table [2] shows the fitted coefficients for the best-fit 
template combination. Fitting the coefficients from each 
template to a power law in frequency yields spectral in- 
dices —2.5 ±0.1 for emission traced by the 408 MHz sur- 
vey and —2.0 ± 0.1 for emission traced by the Cii map. 

The best template model uses the 408 MHz survey to 
trace synchrotron emission and the FIRAS map of Cii 
line emission to trace free-free emission. We have re- 
peated the template fits using different tracers for syn- 
chrotron or free-free emission, obtaining broadly similar 
results. In all cases, the best choice for free-free tem- 
plate is either the Cii map or the thermal dust map. We 
note that spatial structure in the ARCADE 2 maps is 
dominated by the Galactic plane, where extinction cor- 
rections for free-free maps derived from Ha emission are 
worst. The ratio of emission traced by the Cii free- free 
template to emission traced by the 408 MHz synchrotron 
template is (Tff/Tsynch) = 0.16 ±0.01 evaluated for lati- 
tudes l&l < 10° in the lowest ARCADE 2 channel at 3.15 
GHz. The ratio of free-free emission to the total Galactic 
plane emission is (Tfi/Ttotai) = 0.10 ± 0.01. 

4. GALACTIC POLAR CAP TEMPERATURE 

The sky temperature measured by ARCADE 2 in- 
cludes contributions from Galactic, extragalactic, and 
cosmic sources. The 2.7 K cosmic microwave back- 
ground dominates the measured temperatures. Galac- 
tic emission is dominated by synchrotron emission, with 
a smaller contribution from free-free sources. The inte- 
grated contribution of similar synchrotron and free-free 
emission in external galaxies constitutes an isotropic ex- 
tragalactic radio background. Accurate measurement of 
the extragalactic or cosmic backgrounds requires reliable 
determination of both the spatial structure and the zero 
level of the Galactic emission. 

Although the cosmic microwave background may be 
distinguished from Galactic or extragalactic radio emis- 
sion by the different frequency dependences, spectral fit- 
ting alone can not distinguish Galactic emission from 



an extragalactic component of similar spectral behav- 
ior. The template model above reproduces the spatial 
structure of the observed Galactic emission, but is in- 
sensitive to emission described by an additive constant 
in each map. If the total Galactic brightness were de- 
termined along some fiducial line of sight, we could add 
a monopole term to the template model to match the 
total Galactic brightness along that line of sight. The 
adjusted model would then characterize Galactic emis- 
sion along all other lines of sight. 

The north and south polar caps (latitude |6| > 75°) 
provide convenient reference lines of sight. We use two 
independent methods to estimate the total Galactic emis- 
sion toward the Galactic polar caps. The first method 
relies solely on the spatial morphology, assuming a plane- 
parallel model to estimate the total Galactic emission at 
the polar caps. A second, independent method compares 
the spatial morphology of the ARCADE 2 and other ra- 
dio maps to template maps of atomic line emission in 
order to estimate the amount of radio emission per unit 
line emission. Since line emission can unambiguously be 
attributed to the Galaxy, the total Galactic radio emis- 
sion toward the polar caps may then be estimated by 
scaling the total line emission at the caps by the observed 
ratio of radio to line emission. Both methods yield simi- 
lar values for the total Galactic emission in the polar cap 
regions. 



4.1. Polar Cap Temperature From csc(6) 




Cosecant( |b| ) 



Fig. 3. — Temperature of the 408 MHz survey (top) and AR- 
CADE 3.15 GHz channel (bottom) binned by csc(f)) for latitudes 
b > 10° (Northern hemisphere). Sky temperatures include the 
CMB and any extragalactic background. The scatter about the 
best-fit line results from the higher-order spatial structure in the 
maps. For clarity, the plotted uncertainties have been inflated by 
a factor of 20 (408 MHz survey) or 5 (ARCADE); the actual sta- 
tistical uncertainties for the mean of each bin are smaller than the 
symbol size. The data are consistent with a plane-parallel model 
whose measured slope provides an estimate of the Galactic emission 
at the poles. 
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Fig. 4. — Galactic emission TqIu) towards the North galactic 
pole, derived from a cosecant fit to the spatial structure (open 
circles) and the radio/Cll correlation (filled circles). Both methods 
agree and are consistent with a single power law over the frequency 
range 22 MHz to 10 GHz. The fitted Galactic temperature towards 
the polar caps can be combined with the template model of the 
spatial structure to fully specify the Galactic model. 

A commonly used method to estimate total Galactic 
brightness models the Galaxy as a simple plane-parallel 
structure and fits the spatial distribution of the radio 
continuum at Galactic latitude |6| > 10° to the form 



TA{iy,p) = c{iy)+TG{i^)csc{\b\) 



(5) 



where b is the Galactic latitude of each pixel p and 
the constant c{i>) accounts for the contribution from the 
CMB or extragalactic backgrounds at each frequency ly. 
Figure [3] shows typical results for the northern hemi- 
sphere, both for the lowest frequency ARCADE 2 chan- 
nel and the 408 MHz survey. The high-latitude structure 
of the radio sky is consistent with a plane-parallel model, 
whose slope defines the antenna temperature TQiv) of 
Galactic emission toward the Galactic poles. 

We repeat the cosecant fit (Eq. [5]) independently 
for the northern and southern hemispheres, using each 
of the ARCADE 2 sky maps as well as a selec- 
tion of l ower- frequency rad io surveys. Surveys at 
22 M Hz (iRoger et al. I I1999D . 45 MHz (iMaeda et al. " 



1999t lAlvarez et al. I I1997D. 408 MHz (fHaslam et ah 



1981. and 1420 MH z (jReich. Testori. fc Reich 1 120011 : 



Reich &: Reich I I1986D have full or nearly full sky cov- 
erage at frequencies where the sky brightness is domi- 
nated by Galactic radio emission. For each sky map, we 
compute the cosecant slope Taiv) and its uncertainty, 
including both statistical and calibration uncertainties. 
Although the binned data are broadly consistent with a 
plane-parallel model, higher-order structure causes the 
binned data to show greater scatter about the best-fit 
slope than would be expected given the formal statisti- 
cal uncertainty in the mean of each bin. We account for 
higher-order structure by inflating the statistical uncer- 
tainty to force to unity per degree of freedom about 
the best-fit line. We then add the calibration uncertainty 
(of order 10% for the low-frequency radio surveys) in 
quadrature with the (inflated) statistical uncertainty to 
derive the total uncertainty in the cosecant slope for each 
sky survevy. 

Figure 0] shows the Galactic temperature Tg{v) to- 
wards the north polar cap derived from the spatial mor- 



phology of the maps. Although the sky has signifi- 
cant longitudinal substructure within each latitude bin, 
the binned data for each map strongly supports the 
gross morphology dominated by a plane-parallel struc- 
ture. The derived polar cap temperature is consistent 
with a single power law over the frequency range 22 MHz 
to 10 GHz. 

4.2. Radio/Cll Correlation 

Maps of line emission provide a valuable tracer of 
Galactic emission. Line emission from the diffuse inter- 
stellar medium can be measured over the full sky, has 
a well-defined amplitude along every line of sight, and 
does not suffer from extragalactic contamination. Sev- 
eral lines could be used to trace Galactic microwave emis- 
sion. Ha emission from 3-2 transition in neutral atomic 
hydrogen has been mapped over the full sky and is a well- 
established tracer of free-free emission. However, Ha 
emission predominantly traces emission from the warm 
ionized medium, and may not accurately trace emission 
from the brighter synchrotron component. In addition, 
Ha emission suffers from dust extinction in the Galactic 
plane making it less reliable in regions where the Galac- 
tic radio emission is brightest. The 21 cm fine-structure 
line from neutral hydrogen has also been mapped over 
the full sky. However, since HI emission originates from 
the neutral component of the interstellar medium, it is 
unlikely to trace microwave emission from ionized re- 
gions. The Cii line at 158 /im wavelength from singly 
ionized carbon can be used to trace diffuse radio emis- 
sion. It is optically thin even in the plane, does not suffer 
from dust extinction, and is an important cooling mech- 
anism for the diffuse interstellar medium. It has been 
mappe d over the full sky by the COBE /FIRAS instru- 
ment (Fixs en. Bennett, fc Mather I [l999l ). 

Figure [S] shows the correlation between the ARCADE 
3.15 GHz map and the Cii 158 /im map. The Cn map 
has been smoothed to angular resolution 11?6 to match 
the ARCADE resolution. The data segregate into two 
tracks, an upper track from emission near the Galactic 
center and a lower track from emission near the Cygnus 
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Fig. 5. — Correlation between the ARCADE 3.15 GHz inten- 
sity and the Cll 158 ^m atomic line in units nW ra~^ sr~^ . The 
data clearly divide into two regions each with radio intensity pro- 
portional to the square root of the Cll intensity. The upper track 
originates from pixels near the Galactic center, while the lower 
track originates from pixels near the Cygnus region. 
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Fig. 6. — Correlation between the full-sky 408 MHz survey and 
the Cll 158 fim atomic line in units nW sr~^. The bifur- 

cation observed in the limited ARCADE 2 sky coverage persists 
over larger regions of the sky. We fit the data to two lines and 
iteratively assign each point to one line or the other (see text for 
details). Colors indicate the final assignments. 

region. Both tracks show radio emission proportional to 
the square root of the Cii intensity. This would be ex- 
pected if the radio emission were proportional to the den- 
sity n in the interstellar medium (e.g. from synchrotron 
emission), while the coUisionally excited Cii intensity 
were proportional to n^. Microwave free-free emission, 
proportional to n^, would result in a curved track in 
Figure El The absence of any such curvature provides 
additional evidence that free-free emission is faint com- 
pared to synchrotron at the ARCADE 2 bands. Note, 
though, that if a separate template is provided for the 
synchrotron component, the Cii map may be used to 
trace the fainter free-free component (at least at the AR- 
CADE angular resolution). 

The prominent segregation in Fig [5] results at least in 
part from the limited sky coverage. The brightest regions 
in the ARCADE 2 maps occur at the edge of the Galactic 
center region and again at the edge of the Cygnus region. 
The ratio of Cii intensity to the far-infrared continuum is 
known to differ by a factor of two between these regions 
(jFixsen. Bennett. &: Mather1ll999( l. We investigate the 
extent to which the ratio of Cii intensity to radio con- 
tinuum emission varies by performing a similar analysis 
on radio surveys at 22 , 45, 408, and 1420 MHz. AU 
exhibit evidence for similar segregation, although not as 
pronounced as for the smaller ARCADE 2 sky coverage. 

Figure [6] shows the correlation between the Cii map 
and the 408 MHz survey, after smoothing the 408 MHz 
survey with a 5° FWHM Gaussian to approximate the 
7° tophat FIRAS beam. We quantify the tendency for 
the data to segregate into two tracks by fitting two lines 
to the correlation plot. After a first guess at the fit pa- 
rameters (intercept and slope for each line), we assign 
each point to one line or the other based on the dif- 
ference in antenna temperature between that point and 
either line. We then re-compute the fit parameters for 
each line using only those points assigned to that line, 
and then re-assign all points based on the new fit pa- 
rameters. Several iterations suffice to produce a stable 
solution independent of the initial guess. The segrega- 
tion is statistically significant: fitting the 2538 points to 




Fig. 7. — Pixel masks in Galactic coordinates showing the 
location of emission components from Fig. [6] Pixels with higher 
radio/Cll correlation (red) lie near the Galactic center or along the 
North Galactic spur (radio Loop I). Pixels in blue have radio/Cll 
correlation smaller by a factor of roughly 2. Pixels in black have 
poor signal to noise ratio in the Cll map and arc not used. 

two lines reduces the by a factor of 4 compared to a 
single line fit to all points. 

The segregation in the correlation plot can be related 
to known structures in the radio sky. The points in Fig. [S] 
are colored to highlight the segregation into two compo- 
nents. Figure[7]maps each point in Galactic coordinates. 
Points associated with higher radio/Cii ratio lie toward 
the Galactic center and the North Galactic Spur (ra- 
dio Loop I) . Both are regions with enhanced synchrotron 
emission. 

4.3. Polar Cap Temperature From Cll Correlation 

The observed correlation between radio emission and 
Cii intensity allows a well-defined determination of the 
associated Galactic radio emission toward the Galac- 
tic poles. We use the pixel masks defined by the 408 
MHz/Cii correlation (Figure [7]) to fit each radio survey 
to the form 

2 

TA{iy,p) = ^6,(i.) +a,(i/)( Icip) )°■^ (6) 

1=1 

where Ta is the antenna temperature at frequency v, Ic 
is the Cii intensity, and p denotes the pixels within one 
of the two spatial masks. 

The intercepts bi for the two masks include contribu- 
tions from extragalactic sources (including, notably, the 
CMB monopole) and can not uniquely determine abso- 
lute level of Galactic emission. We instead define the 
Galactic emission at the polar caps by extrapolating the 
observed radio/Cii correlation to high latitude: 

TG(i^) = a(z.)4",^p, (7) 

where /cap is the mean intensity of the Cii map at the 
polar caps (Galactic latitude |5| > 75°). 

The two spatial masks separate the sky into regions 
with higher or lower radio emission per unit Cii inten- 
sity. The difference between these two components is less 
clear at high latitude where both the radio emission and 
Cii intensity become weaker. The spatial segregation of 
the two components suggests that the south polar cap re- 
gion is associated with the weaker radio/Cii component 
(Figure [7|) . The brighter radio regions associated with 
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TABLE 3 
Radio/Cii Correlation Slope 



Frequency 


Correlation Slope a (mK [ nW m sr 


.11-0.5 ■) 


(GHz) 


Mask 1 


Mask 2 


Mean 


0.022 


(1.12 ± 0.06) X W 


(0.74 ±0.02) X lO* 


(0.93 ±0.19) X 10* 


0.046 


(2.74 ±0.07) X 10^ 


(1.56 ±0.03) X 10^ 


(2.15 ±0.63) X 103 


0.408 


18.1 ±0.3 


8.5 ± 0.1 


13.3 ± 5.0 


1.420 


(5.11 ± 0.08) X 10-1 


(2.94 ±0.04) X 10-1 


(4.02 ± 1.10) X 10-1 


3.195 


(5.63 ± 0.12) X 10-2 


(2.88 ±0.09) X 10-2 


(4.25 ± 1.38) X 10-2 


3.300 


(5.05 ± 0.11) X 10-2 


(2.61 ±0.08) X 10-2 


(3.82 ± 1.20) X 10-2 


8.15 


(5.73 ± 0.37) X 10-^ 


(1.99 ±0.36) X 10-=* 


(3.86 ± 1.87) X 10-3 


8.33 


(6.32 ± 0.34) X 10-3 


(2.59 ±0.26) X 10-3 


(4.46 ± 1.86) X 10-3 


9.72 


(3.68 ± 0.15) X 10-=* 


(2.37 ±0.16) X 10-3 


(3.03 ±0.66) X 10-3 


10.15 


(3.67 ±0.16) X 10-=* 


(2.17 ±0.18) X 10-3 


(2.92 ±0.75) X 10-3 



the North Galactic spur, however, extend to high lati- 
tude so that the mean radio brightness associated with 
Cii emission toward the north polar cap region proba- 
bly lies between the limiting cases established by the two 
masks. We therefore use the mean of the slopes from 
the two masks to estimate the Galactic radio emission 
at the north polar cap, with uncertainty broad enough 
to bracket the two cases. The mean value thus derived 
is consistent with similar analysis fitting a single slope 
to the full sky. Within the restricted ARCADE 2 sky 
coverage, however, the value of a single slope fitted to 
the entire observed region could vary significantly as the 
sky coverage shifts slightly to include more or less of the 
Galactic center region. The mean of the two observed 
slopes is less susceptible to variations in the ARCADE 2 
sky coverage. 

Table[3]shows the correlation slopes ai{v) fitted to each 
mask for both the ARCADE 2 data and lower-frequency 
sky surveys. Figure |4] compares the Galactic temperature 
towards the north polar cap, derived from mean Cii cor- 
relation, to the value derived from the cosecant fit. Both 
methods provide similar estimates for the Galactic com- 
ponent of the total brightness temperature towards the 
north polar cap. This Galactic component is consistent 
with a single power law 

TG{v)^TG.x{y/yoY (8) 

with spectral index (3 = —2.55 ± 0.03 and amplitude 
Tcai = 0.498 ± 0.028 K at reference frequency vq ^ I 
GHz. 

The analysis above uses a single map (Cii line emis- 
sion) to trace Galactic emission from the diffuse inter- 
stellar medium. If a significant fraction of Galactic radio 
emission originated from a component of the interstellar 
medium not well sampled by Cii emission, the result- 
ing correlation would under-estimate the total Galac- 
tic emission toward the polar cap. We test for addi- 
tional microwave emission from other components of the 
interstellar medium by repeating the correlation analy- 
sis using a simultaneous fit to three line maps chosen 
to sample different components of the diffuse interstel- 
lar medium: the (square root of the) COBE/ FIRAS 
map of Cii emission, a map of Ha emission (FinkbeinerJ 
l2003fl. and the Leiden/Arg entine/Bonn Galactic HI Sur- 
vey ()Kalberla et al. II2005I ). The simultaneous fit pro- 
duces results nearly identical to a fit using just the Cii 
map, shifting the derived Galactic polar cap temperature 
by 11 ± 21 mK at reference frequency 1 GHz. Since the 
limited ARCADE 2 sky coverage increases the effect of 



covariance between the different line maps, we use the 
results derived from the single Cii correlation. 

4.4. Composite Galactic Model 

Removing Galactic emission from the ARCADE 2 data 
requires a model for both the spatial structure and the 
total Galactic emission. At each ARCADE frequency, we 
generate a full-sky model of the spatial structure using 
full-sky template maps niultipled by the coefficients fitted 
within the region observed by ARCADE 2 (Tabled. We 
then add a constant to the template model to match the 
total Galactic emission derived from the Cii and csc(|&|) 
fits towards selected reference positions. 

We independently model the total Galactic emission 
toward the north and south Galactic polar caps > 
75°). For each reference position, we compute the to- 
tal Galactic emission associated with either the plane- 
parallel structure (Eq. [5]) or the Cii emission (Eq. [T]). 
We apply each method to the ARCADE 2 and radio sur- 
veys, and parameterize the resulting multi-frequency re- 
sults using a power-law model (Eq. [8]) . As a cross-check, 
we also estimate the total Galactic emission for the cold- 
est patch in the Northern hemisphere, consisting of all 
pixels within 15° of 6 = 48°, I = 196° (north of the Galac- 
tic anti-center). Since by definition a mid-latitude cold 
spot is inconsistent with a plane-parallel structure, the 
estimate for this position is based only on the radio/Cii 
correlation. Table [4] shows the estimated Galactic nor- 
malization Tcai and spectral index f} for each reference 
position. We obtain similar results using either the spa- 
tial morphology or the radio/Cii correlation. 

We fix the offset of the template model (Eq. [3]) by 
computing the temperature of the template model to- 
wards these same regions, then adding a constant ao{y) 
to the template model at each frequency v to force the 
template model to match the power-law model along 
these lines of sight. The resulting composite model 
of Galactic emission can be used to correct the cali- 
brated time-ordered data in order to estimate the CMB 
monopole tem perature and t he extragalactic radio back- 
ground (Fixse n et al. |[2b08f ). The offsets derived from 
the three independent regions agree within 5 mK at 3 
GHz, which we adopt as the uncertainty in the zero level 
of composite Galactic model. 

The spectral index derived for the Galactic emission to- 
ward the polar caps or coldest region, (3 k, —2.57 ± 0.03, 
is consistent with the spectral index (3 — —2.5 ± 0.1 de- 
rived from the synchrotron template fit only to ARCADE 
2 data, indicating that high-latitude Galactic emission is 
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TABLE 4 

Galactic Emission" Along Selected Lines of Sight 



Parameter 


I'ochniquo 


North Polar Cap 


South Polar Cap 


Coldest Patch 


^Gal (K) 


Hadio/Cll 
csc(|6|) 

Weighted Mean 


0.492 ± 0.095 
0.499 ± 0.030 
0.498 ± 0.028 


0.303 ±0.054 
0.366 ±0.028 
0.353 ±0.025 


0.188 ± 0.130 
0.188 ±0.130 


/3 


Radio/Cll 
csc(|6|) 

Weighted Mean 


-2.53 ±0.07 
-2.56 ±0.04 
-2.55 ±0.03 


-2.59 ± 0.06 
-2.65 ± 0.05 
-2.63 ± 0.04 


-2.56 ±0.12 
-2.56 ±0.12 



" Galactic emission Tg{u) = Tqi^i{u/uo)^ with reference frequency uq = 1 GHz. 



dominated by synchrotron emission at frequencies below 
10 GHz. A synchrotron spectral index of -2.57 is sig- 
nificantly flatter than the values derived from measure- 
ments at higher frequencies, but is consistent with the 
majority of measurements below 10 GHz (see, e.g., the 
tabulation in Rogers fc Bowman . (,,2008.)). Notable ex - 
cept ions are the measureni ents by iTartari et al. I (|2008D 
and iPlatania et al. I (|19980 . which prefer /3 « -2.8 for 
the synchrotron component. Both of these results are 
derived from a single strip at declination S — +42° and 
may not be representative of larger regions of the sky. 

5. SPINNING DUST 

The ARCADE 2 data can also be used to search 
for additional microwave emission associated with in- 
terstellar dust. The microwave sky is known to 
contain a component spatially correlated with far- 
infrared dust emission but not with synchrotron- 
domin ate d surveys at 408 or 1420 MHz ( Koeut et al 
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1996a|lbl: Ide Oliveira-Costa et al. I Il997t 
1993 iBennett et al. Il2003[ ). The freauencv spectrum of 



the correlated component differs markedly from the emis- 
sion spectrum of thermal dust, showing a spectral index 
f3 w —2.2 in antenna temperature from 20 to 50 GHz. 

Two main candidates have emerged to explain this 
"anomalous" correlated emission component. Electric 
dipole emission from a population of small, rapidly 
rotating dust grains will produce a spectrum with 
a broad peak in t he fr equ ency range 20-40 GHz 
Draine fc Lazarianl Il998f ). iMiville-Deschenes et al. I 
2008f ) analyze data from the Wilkinson Microwave 



Anisotropy Probe (WMAP) mission and conclude that 
spinning dust accounts for the majority of emission in 
the Galactic plane at 22 GHz. An alternative model 
explains the correlated emission as flat-spectrum syn- 
chrotron emission associated with star fo rmation activity 
(|Bennett et al. Il2003t [Gold et al. II2008D . Energy losses 
as electrons propagate from their origin steepen the syn- 
chrotron index away from these sources. Since a flat- 
spectrum component will increasingly dominate at higher 
frequencies, the spatial morphology of synchrotron emis- 
sion at higher frequencies should increasingly resemble 
thermal dust emission from the same star formation ac- 
tivity. 

Efforts to distinguish between these models have 
largely focused on the frequency spectrum of the dust- 
correlated component inferred by correlating a map trac- 
ing thermal dust emission against microwave maps at dif- 
ferent observing frequencies. The correlation with ther- 
mal dust emission is seen to peak at frequencies between 
20 and 30 GHz, falling at frequencies below 22 GHz 
(|de Oliveira-Costa et al. Il2004l : [Fternandez-Cerezo et al. I 



120061 : iHildebrandt et al. l[200l . This falling spectrum 
at low frequencies has been cited as evidence favoring 
spinning dust models. 

The decrease in dust-correlated emission at frequencies 
below 20 GHz, although consistent with spinning dust 
emission, does not necessarily support such models over 
the flat-spectrum synchrotron alternative. Both models 
predict weaker correlation with thermal dust emission at 
frequencies below 20 GHz. For the spinning dust model, 
the weaker correlation results from the lower amplitude 
of the spinning dust emission. The flat-spectrum syn- 
chrotron model, however, also predicts a weaker correla- 
tion at lower frequencies, since by construction the spa- 
tial morphology of the sky in this model varies smoothly 
from the WMAP data at 22 GHz (observed to correlate 
well with the thermal dust morphology) to the Haslam 
map at 408 MHz (observed to correlate only weakly with 
thermal dust emission). Model discrimination based on 
cross-correlations with thermal dust emission must rely 
on rigorous statistical comparison of competing models, 
and not simply on the general trend toward weaker cor- 
relations at lower frequencies. 

A simpler test relies on the absolute spectrum of Galac- 
tic emission. At frequencies below 30 GHz, thermal dust 
emission is negligible, so that the total Galactic emis- 
sion becomes a superposition of free-free, synchrotron. 




Frequeiicj' (GHz) 

Fig. 8. — Synchrotron and spinning dust emission for two models 
with different spinning dust amplitude. The red (blue) curves show 
emission from a model with high (low) dust normalization. Dashed 
lines show the dust emission, dotted lines show synchrotron emis- 
sion, and solid lines show the combined emission from each model. 
The amplitude of the combined emission is fixed at both 22 GHz 
and 408 MHz. Measurements below 10 GHz are dominated by 
synchrotron, which can be used to infer the dust normalization: 
models with more dust emission at 22 GHz have lower synchrotron 
emission. Grey bars indicate the ARCADE 2 frequency bands. 
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and spinning dust emission. Both free-free and syn- 
chrotron emission increase monotonically at lower fre- 
quencies. Spinning dust, in contrast, decreases in ampli- 
tude below 20 GHz, so that the spectrum of the com- 
bined Galactic emission below 20 GHz can be used to 
place limits on the contribution of spinning dust without 
reference to detailed spatial correlations with a thermal 
dust template. 

Figure [5] illustrates the concept. Consider Galac- 
tic emission consisting of a superposition of syn- 
chrotron and spinning dust (ignoring for the moment the 
smaller contribution from free- free emission). Following 
iMiville-Deschenes et al. I (|2008l l , we specify the spinning 
dust amplitude as a fraction of the total Galactic bright- 
ness at 22 GHz, assumed here to represent the peak in the 
dust spectrum. As the dust normalization is increased, 
the synchrotron amplitude at 22 GHz must decrease to 
keep the total emission constant. The spinning dust spec- 
trum falls rapidly, so that emission below 10 GHz is dom- 
inated by the synchrotron component. A model with 
high dust normalization (red curves in Fig. [5]) will thus 
have fainter synchrotron emission at frequencies of a few 
GHz, while a model with less spinning dust (blue curves) 
will have brighter synchrotron emission. This shift in the 
amplitude of Galactic emission at frequencies below 10 
GHz is a sensitive test for spinning dust, without resort 
to detailed spatial correlations. 

We implement this test using a simple model of Galac- 
tic emission. We assume that the spatial distribution 
of spinning dust emission is traced by a template map 
of thermal dust emission, and fix the amplitude of as- 
sociated spinning dust emission by scaling the template 
map so that the re-scaled dust map forms a specified 
fraction of the Galactic plane intensity measured by the 
WMAP 22 GHz map. This allows a simple normalization 
of spinning dust in terms of its relative contribution to 
the total Galactic plane intensity at 22 GHz. We use the 
COBE/DIRBE 240 ^m map of thermal dust emission 
(|Reach et al. Ill996[ ) as the thermal dust template since 
it is dominated by thermal dust emission but unaffected 
by extinction in the plane. For specificity, we define the 
Galactic plane mask using all pixels lying within the AR- 
CADE 2 observation pattern with latitude |6| < 20°, and 
use this mask for all computations so that the model re- 
sults are not affected by differing sky coverage between 
ARCADE 2 and other surveys. 

We then model the Galactic emission spectrum as 
follows. We first correct the WMAP 22 GHz map 
and the Haslam 408 MHz map by subtracting the 
WMAP max imum-entropy model of free-free emission 
(|Gold et al. 1 120081 using a spectral index /3ff = —2.15. 
The corrected maps then contain only synchrotron and 
(possibly) spinning dust emission. We assume that the 
spinning dust spectrum peaks at 22 GHz and scale the 
nor malized spinning dust map to lower frequencies using 
the iDraine fc Lazarian I (|1998[ ) model for the warm neu- 
tral medium. After correction for free-free and spinning 
dust emission, the 22 GHz and 408 MHz maps contain 
only synchrotron emission, which we use to define the 
synchrotron amplitude and spectral index for each pixel. 

The resulting model (synchrotron, free-free, and spin- 
ning dust) can be used to estimate the combined Galac- 
tic emission at frequencies between 408 MHz and 22 
GHz. We compare the model to the ARCADE 2 data 



by smoothing the model map of combined emission to 
the ARCADE 2 angular resolution, and then computing 
the mean Galactic emission of the smoothed model for 
all pixels within the ARCADE 2 galactic plane mask. 
Figure [9] shows the result. If spinning dust is negligible, 
the synchrotron contribution at 22 GHz is maximal and 
the model approximates the flat-spectrum synchrotron 
model. As the spinning dust amplitude increases, the 
synchrotron contribution at 22 GHz decreases and the 
synchrotron spectral index steepens. The combination 
of lower synchrotron amplitude and falling dust spec- 
trum combine to lower the total model emission across 
the ARCADE 2 frequency bands. The ARCADE 2 data 
lie below the model prediction for no spinning dust, and 
are consistent with spinning dust contributing 0.4 ± 0.1 
of the total K-band Galactic plane emission. 

WMAP is a differential instrument and is insensitive 
to any monopole emission component. The zero level of 
the 22 GHz map is s et using the cosecant dependence 
on Galactic latitude fflinshaw et al. Il2008f l. The 408 
MHz survey and the ARCADE 2 data, however, both 
include monopole contributions. To prevent discrepant 
treatment of the map zero levels from affecting the model 
predictions, we remove a monopole from the 408 MHz 
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Fig. 9. — Mean intensity of Galactic plane emission for the 
ARCADE 2 data, compared to model predictions with and without 
spinning dust (see text). Top panel: Model predictions and data. 
The upper solid curve shows the model prediction for no spinning 
dust, averaged over the ARCADE 2 sky coverage with < 20°. 
The lower solid curve shows the model prediction for spinning dust 
amplitude equal to 60% of the total Galactic plane emission at 22 
GHz. The dotted curve shows the best fit to the ARCADE 2 data. 
The spectrum of emission from spinning dust is assumed to peak 
at 22 GHz, so a higher spinning dust fraction at 22 GHz produces 
lower total emission at the ARCADE 2 frequencies. The bottom 
panel shows the same data, normalized by dividing each point by 
the model prediction for no spinning dust. Galactic plane emission 
observed by ARCADE 2 is consistently fainter than expected for 
a model with no spinning dust, and is consistent with spinning 
dust contributing 0.4 it 0.1 of the total Galactic plane emission at 
reference frequency 22 GHz. 
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survey and the ARCADE 2 sky maps, and re-set the 
zero level of each map using a csc(|&|) fit to the spa- 
tial structure in each map to match the processing of 
the WMAP data. We restrict analysis to pixels at low 
Galactic latitude (|6| < 20°) where Galactic emission is 
brightest so that uncertainties in the zero level have min- 
imal effect. The largest uncertainty in the analysis is 
the correction for free-free emission. We use the WMAP 
maximum-entropy model of free-free emission, derived 
assuming that spinning dust contributes negligibly to the 
total Galactic emission. We test the sensitivity of the re- 
sult to the free-free normalization by repeating the anal- 
ysis using the same maximum-entropy model scaled by 
a constant normalization factor. Changing the free-free 
normalization by a factor of 2 affects the best-fit spin- 
ning dust normalization by approximately 0.08, and is 
included in the total quoted uncertainty. 

6. CONCLUSIONS 

We use the ARCADE 2 absolutely calibrated obser- 
vations of the sky to model Galactic emission at fre- 
quencies 3, 8, and 10 GHz and angular resolution 11?6. 
TT plots of the binned sky maps show net spectral in- 
dex P = -2.43 ± 0.03 between 3.3 and 8.3 GHz and 
f3 = -2.47 ± 0.02 between 3.3 and 10.3 GHz, consistent 
with a superposition of synchrotron and free-free emis- 
sion. The spatial structure in the maps can be described 
using two spatial templates, the Haslam 408 MHz sur- 
vey to trace synchrotron emission and the COBE/FIRAS 
map of Cii emission to trace free-free emission. Fitting 
these templates to the ARCADE 2 maps yields spec- 
tral index /3synch = —2.5 ± 0.1 for emission traced by 
the synchrotron template, with free-free normalization 
(Tff/Ttotai) = O.lOiO.Ol evaluated for latitudes |&| < 10° 
in the lowest ARCADE 2 channel at 3.15 GHz. 

The template model only specifies Galactic emission 
up to an additive constant. We fully specify the Galac- 
tic model by computing the temperature of the template 
model toward three reference locations (north and south 
polar caps plus the coldest patch in the Northern hemi- 
sphere) , and adding a constant to the template model to 
match the Galactic temperature derived by other means. 
One estimate of the Galactic temperature models the 
Galaxy as a simple plane-parallel structure and derives 
the polar cap temperatures from the slope of the Galac- 
tic emission binned by the cosecant of Galactic latitude. 
A second, independent method computes the correlation 
between radio emission and Cii line emission, then esti- 
mates the radio brightness at each reference position by 
multiplying the Cii intensity at that position by the mea- 
sured radio/Cii correlation slope. Both methods yield 
similar estimates for the total Galactic temperature to- 



wards each reference location. We extend the analysis to 
full-sky surveys at lower frequencies and find that Galac- 
tic polar cap temperature from both methods is consis- 
tent with a single power law over the frequency range 22 
MHz to 10 GHz, with spectral index /? = -2.55 ± 0.03 
and normalization 0.498 ±0.028 K at reference frequency 
vo = 1 GHz. 

ARCADE2 produces well-calibrated maps of Galactic 
emission. We use these maps to test for contributions 
from spinning dust near the Galactic plane. Previous 
tests for spinning dust at frequencies below 20 GHz have 
used only the spatial correlation with template maps 
tracing thermal dust emission. The frequency depen- 
dence of the dust-correlated component is not a strin- 
gent test for spinning dust, since a weaker correlation 
at lower frequencies is expected both for models with 
significant spinning dust contributions as well as mod- 
els with flat-spectrum synchrotron but no spinning dust. 
We use a simple model of the total Galactic emission to 
predict the mean intensity of the Galactic plane in the re- 
gion observed by ARCADE 2. The model normalizes the 
spinning dust contribution at 22 GHz, and computes the 
expected total contribution from free-free, synchrotron, 
and spinning dust as a function of frequency and spinning 
dust normalization. The ARCADE 2 data consistently 
show less emission in the Galactic plane than a model 
with no spinning dust, and are consistent with spinning 
dust contributing 0.4±0.1 of the Galactic plane emission 
at 22 GHz. 
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